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Abstract 



As a first step towards constructing chiral models on the lattice with staggered fermions, 
we study a U(l) model with axial-vector coupling to an external gauge field in two di- 
mensions. In our approach gauge invariance is broken, but it is restored in the classical 
continuum limit. We find that the continuum divergence relations for the vector and axial- 
vector currents are reproduced, up to contact terms, which we determine analytically. The 
current divergence relations are also studied numerically for smooth external gauge fields 
with topological charge zero. We furthermore investigate the effect of fluctuating gauge 
transformations and of gauge configurations with non-trivial topological charge. 
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1 Introduction 



A consistent non-perturbative formulation of a chiral gauge theory is much desired. Within 
perturbation theory the quantization of chiral gauge theories has been studied and it is 
claimed to be solved satisfactorily |]J , although a gauge invariant regulator does not appear 
to exist. Various proposals have been made for a non-perturbative formulation of chiral gauge 
theories on the lattice. A problem is the doubling phenomenon: a gauge invariant model on 
a regular lattice is necessarily non-chiral in the sense that each fermion is accompanied by 
extra degrees of freedom, the so-called species doublers, which couple with opposite chiral 
charge to the gauge fields and render the theory vector-like. Most of the currently existing 
lattice proposals try to eliminate the unwanted species doublers either by making them very 
heavy or by tuning their interactions to zero (for an overview see ref. 0). One can also use 
the doublers as physical degrees of freedom, with the staggered fermion method [3-5], which 
is the strategy followed in this paper. 

As in perturbation theory in the continuum, the regulated lattice model violates gauge 
invariance, but it is restored in the classical continuum limit. The question is how to restore 
it in the quantum theory. One can mimic the continuum methods closely, by attempting 
non-perturbative gauge fixing and adding counterterms to restore gauge invariance |J. An 
alternative approach focuses on a possible dynamical restoration of gauge invariance |3|, 
The question of dynamical gauge symmetry restoration will be investigated in an other 
publication |7] . In this paper we shall test the staggered method in a two-dimensional axial- 
vector model, where the staggered fermion fields are coupled to external gauge fields. A 
preliminary account of this work has already been presented in ref. ||. 



2 The target model 

Our continuum target model is given by the following euclidean action in two dimensions 

S = - J d 2 x {jfa^dp + i^A^ + rm/nf)} , (2.1) 

where is an external gauge field and 75 = — 27172- For the purpose of numerical simu- 
lations and also for the use in tests we have added a bare mass term for the fermions with 
mass parameter m. In two dimensions the axial- vector model ( |2.1| ) can be rewritten in vector 
form with a Majorana mass term, e.g. by a charge conjugation transformation on the right 
handed fermion fields. So it is not really a chiral gauge theory, and for m = it is equivalent 
to the Schwinger model. However, the technical aspects of '75' in its lattice version are very 
similar to truly chiral gauge theories. The reason for choosing this target model and not 
e. g. a left-handed model, is that the staggered fermion version of the axial model has a 
larger lattice symmetry group. 

For m = the action fl2.1| ) is invariant under the local gauge transformations A^x) — > 
A^x) + i£l(x)d^*(x) } ^{x) -> (n(x)P L + tt*(x)P R )i>(x),^(x) ^^(x)(n(x)P L + n*(x)P R ), 
with Pl,r = (1 T 7s)/2 and Q(x) G U(l). The action is furthermore invariant under the 
global vector symmetry ip(x) — > Qip(x), ip{x) — > tp(x)Q*. The vector and axial-vector 
currents, = iip^^p and = it/)') ^75^ ', satisfy the classical divergence equations 

d^l = 0, (2.2) 
d^tf = 2mJ p , (2.3) 

where J p = iip'jsip is the pseudoscalar density. 



Figure 1: Feynman diagrams for T^ u A (p) , T^ p {p), T AA {p), andT AP (p). 



As is well known these classical relations ( |2.2| ), ( |2.3| ) may be invalidated by quantum 
effects. Requiring gauge invariance for m = 0, eq. ( |2.3|) has to remain valid but ( |2.2| ) 
becomes anomalous and takes the form 

8, Jj = 2zg , g = i-F 12 , (2.4) 

where F 12 (x) = diA 2 (x) — d 2 A 1 (x) is the field strength and q(x) is the topological charge 
density. The topological charge Q is defined by Q = J d 2 x q(x). Our aim is to study the 
divergence equations ( |2.2| ) and ( |2.4j ) in the lattice version of the model. 

Let us briefly review the anomaly structure of the current divergences in the quantum 
theory. We evaluate (d^J^ ,A ) and (J p ) in perturbation theory, by expansion in A^, which 
leads to a series of diagrams in which external A^-lines are attached to a fermion loop. We 
concentrate on the diagrams shown in fig. [IJa, since diagrams with one external line vanish 
and diagrams with more than two external lines are convergent by power counting. We 
shall evaluate these diagrams using a spherical cut-off in momentum space, because this 
regularization has analogies to our staggered fermion formulation introduced later. Such a 
cut-off violates gauge invariance, but the desired result is easily obtained by adding suitable 
contact terms. For a discussion (including e.g. the gauge invariant point splitting method, 
which gives rise to the additional diagrams in fig. |l|b) see ref. ||. 

Although the diagrams in fig. [l|a seem to be logarithmically divergent, the momentum 
cut-off gives a finite answer when it is removed. Fig. la leads to 



T* A (p) = jd 4 xexp(-zxp)(J^(x)J A (0)) A=0 



2tt Jo m z + 2(1 — z)p z ^ 

where C^ A is a contact term to be determined shortly. From this expression we get the two 
Ward identities 



W^u A (p) = ^uaPa + ip^C™ 



1 

2tt' 

-iPuT^ A (p) = 2mT v /{p) - ^e, aPa - ip v C^ A , (2.6) 



where the amplitude 



results from the evaluation of the the J^-J p correlation function at = (see fig. |]a). 
These relations show anomalous terms in both the vector and axial-vector Ward identities. 
By choosing the contact terms CY^ = ±(i/2n)e^ u one may shift the anomaly either to the 
vector (+) or to the axial- vector Ward identity (— ). Since we insist on axial gauge invariance 
the contact term is determined as 

= (2.8) 



giving 



W^ v A {p) = l -t vaVa , -ip u Tj v A (p) = 2mTj P ip) ■ (2.9) 



For fig. la gives 



)A=o 



1 I ' dz (S„y - - z) + m V + 1, + . (2 . 10) 

7T JO 771T + -2(1 — ZjJT 27T 



With the choice 



2vr 

we get from this relation the desired Ward identity 



C 4 = I 2 ' 11 ) 



- ™f W , r-(p) = -± /; ^ m2 + ^"_ z)p2 , (2.12) 



with T AP {p) the J A -J P correlation function at = 0. From the above results for the 
current correlation functions one can derive the divergence equations ( |2.3| ) and ( |2.4| ). 

We shall show in the following that the symmetry properties of the model, as they 
manifest themselves in the current divergence relations ( |2.3| ) and (2^1), can be recovered in 
a staggered fermion version of the model on the lattice. 



3 The staggered fermion model 

We generalize our target model to two flavors, which makes it somewhat simpler to describe 
in the staggered fermion formalism. We introduce 2x2 matrix fermion fields \1/" K and ^> x 
on a two dimensional square lattice, where a and k are Dirac and flavor indices, respectively. 
Using these matrix fields we find after the naive lattice transcription of the two flavor version 
of the target model in eq. (|2.1| ) the following action, 

s = - £ l Tr {^(tf x p L + u p p R )^ x+fl - u x+{iltx {u^p L + u p :p r )^ x } 

XTr{^}, (3.1) 



XfJ, 

m 



with a the lattice spacing. We shall use lattice units, a — 1. The action above would be 
gauge invariant if ^! X K and ^ x would be independent degrees of freedom, and this would 
lead to fermion doublers with opposite chiralities. 

The fermion doublers are situated at the boundary of the Brillouin zone in momentum 
space, i.e. at momenta p^ = ±7r. Consider restricting the momenta of the matrix fields 
such that — 7r/2 < p^ < +tt/2. Then we loose the fermion doublers of the matrix fields, but 
also gauge invariance. It is clear, however, that in the classical continuum limit, where the 
field momenta go to zero (in lattice units), eq. Q3.1| ) goes over into the eq. ( |2.1| ) with two 
flavors. Hence gauge invariance gets restored in this limit. One might think that the cut-off 
in momentum space has to result in a non-local action. However, it is possible to express 
the action in a form that is local, using staggered fermions. 

Staggered fermion fields on the lattice, denoted by the one-component fields Xx and 
X x , do not carry explicit flavor and Dirac labels. These labels are supplied through the 
doubling phenomenon. In the classical continuum limit one recovers the usual Dirac and 
flavor structure. We make the connection with the ^ X K and fy x a by writing p, [5] 

*» = ITS E T* +6 X,+6 , = 173 E(7* +6 fe +b , (3.3) 

b b 

where j x = 7i 1 72 2 and the sum runs over the corners of an elementary lattice square, 
bfj, — 0, 1. In momentum space we have the relation Q 

^ a .(p) = z(p)Y / T aK ,b(p)x(P + ^b), -n/2<p<n/2, (3.4) 

b 

where T aK ^ = X) c ex P(^ C7r )7aK/4 is a unitary matrix and z(p) is a non- vanishing function 
in the restricted momentum interval. This shows clearly that in this restricted interval the 
Fourier components of the matrix fields are independent. 

Having expressed the matrix fermion fields in terms of the independent Xx and x x fields, 
substitution into the action fl3.1| ) leads to a local action. The indices a and Konf and \I/ 
act like Dirac and flavor indices and one can construct staggered fermion models involving 
arbitrary spin-flavor couplings to other fields in a straightforward manner such that the target 
models are recovered in the classical continuum limit. For the Standard Model and Grand 
Unified Theories like SO(10) and SU(5) this can be done such that the staggered fermion 
symmetry group is preserved ||. This invariance is important for reducing the number of 
counterterms needed to get a satisfactory continuum limit [|Hj]. This strategy of coupling 
the staggered fermion spin-flavors has recently been successfully applied to a fermion-Higgs 



model 11 



By working out the trace in ( |3.1| ) one obtains the action in terms of the staggered fermion 
fields, 



$ — ~ 2 X | C V X Z X Vnx+biXx+bXx+b+p, — Xx+b+fiXx+b) 

Xfl { b 



— S{ix\ X/ VlZx+ciVux+cXx+bXx+c+u — Vux+bXx+h+flXx+c) f ~ m X/ XxXx , (3.5) 
b+c=n ) x 

with the abbreviations c^ x = ReU^, s^ x = lmU^ x and n = (1, 1). The sign factors rji x = 1 
and r]2x = (— l) Xl represent the Dirac matrices 71 and 72 and rji2 X = = ( — -O^ 1 

the «75 = 7x72- In the classical continuum limit this action describes two flavors of axially 
coupled Dirac fermions. 



For a one flavor staggered fermion model we would have used only one Grassmann variable 
per site (so-called 'reduced' or 'real' staggered fermions). One then defines the \ x fields 
usually on the even, and the Xx fields on the odd lattice sites. This would require a one link 
mass term instead of the simple one site mass term in ( |3.5| ). The continuum interpretation 
is in this case somewhat more involved |4], [T^| . 

The couplings in the s^ x term in the action ( |3.5| ) are not confined within a plaquette. 
For example, for c = and b = n we have three link couplings. The action is of course not 
unique. According to standard staggered fermion properties, the three link couplings may 
be replaced by one link couplings by shifting the x or X field over two lattice spacings in 
the same direction (even shifts). Such actions are equivalent in the sense that they lead to 
the same classical continuum limit, and in the quantum theory they are expected to be in 
the same universality class. It is instructive to give here a particularly simple alternative to 

S = ~ {^/J-xV^x^iXxXx+p, ~ Xx+jiXx) 

Xfl 

Snx^fivVvx 2 {XxXx+0 ~l" Xz+pXai)} ^ ] m XxXx > (3-6) 



Cfix 



. c /ix-b i S^x — ^ S^x-b i (3-7) 
b b 



obtained by using c = n — b, the identities rj^x+nVux+n — £ /j,uVvx-, n + fi = v + 2/2, and the 
equivalence of e.g. XxXx+n+u+2fi-2b with XxXx+n+o which differ by even shifts. A further 
reduction could be achieved by the replacements c MX — > 1 and — ► s M in the above 
expression, as the resulting model still has the same classical continuum limit. 

A model with all the couplings confined within a plaquette may be called a canonical 



model, as it allows for a canonical construction of the transfer operator |T2|, [L3|]. In the 
following we shall use, however, the original version as written in eq. ( |3.5|) . 

For m = the action (|3.1| ) appears to be invariant under the local gauge transformations, 

Ufix > ^xUfix^x + fl 5 (3-8) 

^ - (si x p l + nzp R )* a , * a -> v x (n x p L + n* x p R ) . (3.9) 

However, this gauge invariance is broken because of the momentum space cut-off on the 
matrix fields. The four components of the ty x and ty x matrix fields are not independent, as 
is evident from their expression in terms of the independent Xx and x x fields. Therefore we 
cannot translate the gauge transformations on \l> and \& to local transformations on x an d 
X, and the action (|3.5| ) lacks gauge invariance. 
The global vector U(l) transformation 

Xx -> exp(iuj)xx , Xx -> Xx exp(-iu) (3.10) 



is an exact symmetry of the actions ( ft.5|) and (3^). For m = there is furthermore a second 



exact global U(l) invariance, the 'U(l) e ' invariance Xx ex 9{^ uje x)Xxi Xx ~^ Xx ex P{i UJ£ x), 
e x = (— l) Xl+X2 , which corresponds to a flavor non-singlet chiral transformation ^> x — > 
cosu ^ x + i sinu; 75^/^75, ^x — * cos to ^ x + i sin to 7s\l/x75- These U(l)xU(l) transformations 
are part of the U(2)xU(2) global invariance of the classical two flavor continuum action. In 
a one flavor staggered fermion model U(1) E would be the only global U(l) symmetry. 

We are not concerned here with the full aspects of flavor symmetry restoration, but 
concentrate on the flavor singlet vector and axial- vector currents, J^ x and J^ x . The latter is 
the gauge current which should be exactly conserved. The former is the U(l) vector current 
which should have the anomaly discussed in the previous section. 



4 Divergence equations on the lattice 



In this section we determine the contact terms in the divergence equations for the gauge 



current J A and a flavor singlet U(l) current 3 in our lattice model ( |3.5|) 



To introduce these currents we generalize ( |3.1| ), ( |3.2| ) to include also an external vector 
field V^x, writing 

IjR = e +iA^-iV^ _ ^ 



U, 



fix i fix 

Since gauge invariance is broken we expect to have to add counterterms to the action to 
restore it in the continuum limit. From the sect. 2 we expect these to have the form 



Set — {\C^A^ X A UX + C^ A V^ X A UX + \C^V^ X V VX ^ . 



(4.2) 



The coefficients C AA and C VA have to be determined such that in the scaling region the 
effective action obtained by integrating out the fermion fields is invariant for m = under 
the gauge transformations ( |3.8|) on A^ x . We have included a C vv term which may be 
determined such that under gauge transformations on the effective action in the scaling 
region suffers only the anomaly and no further symmetry breaking. There is no reason, of 
course, for the numerical values of C AA and C VA to be the same as in sect. 2, where we used 
the spherical cut-off as a regulator. 

The above counterterms may be extended to periodic functions in and V^. E. g. 
writing C AA = t5^ u the C AA term may be replaced by the standard lattice form for a 
gauge field mass term rEj.„(l — cosA^ x ). This replacement could help to reduce the scaling 
violations. In this paper we shall however stay with the form ( [4.2|) . 



The currents Jj ,A are identified from the total action S + S ct by letting A^ — ► A^ + 5 A 



V/j, ~~ > Vfi + SVfi, an d collecting terms linear in 5A^ and 5V^. The current correlation functions 
are obtained by differentiating the effective action with respect to A^ and V^. As we shall 
study (J^ ,A ) only for zero (but for arbitrary A^), we give the currents here for = 0, 



V = -Tr 

° fix 2 

J A = — Tr 

fix 2 



^{U^Pl + U; x P R )V x+fl + ^ X+ ^{U* P L + U, X P R )V X + Cl A A vx , (4.3) 



v.Mu^Pl - u; x p R )* x+fi + * x+(l ^(u; x p L - u, x p r ^ x \ + c aa a vx . 



(4.4) 



where U^ x = U^ x = exp(— iA^). A natural choice for the pseudoscalar density is given by 



J: = tTr 



(4.5) 



The currents in terms of the staggered field x are obtained by inserting the relations ( |3.3D 
into (0)-(0 



J, 



v 

fix 



J, 



fix 



c fix Vfix+biXx+bXx+b+fi + Xx+b+ fiXx+b) 



— Sfix V^x+ciVfix+cXx+bXx+c+fi + Vfix+bXx+b+fiXx+c) 

b+c=n 



Cfix Vl2x+c{?]fix+cXx+bXx+c+fi V/ix+bXx+b+fiXx+c) 

b+c=n 



+ C VA A VX , (4.6) 



+bXx+b+fi — Xx+b+[iXx+bj 



+ C^ A Aux , (4.7) 

b 

J x = 1 J2 Vl2x+bXx+bXx+c ■ (4.8) 
b+c=n 



The above vector current J^ x is not directly related to the exact global U(l) invariance fl3.1U| ) 



which we mentioned in the previous section, because the prescription ( |4.1|) does not make 
this symmetry an exact local symmetry. 

To obtain the divergence relation of the conserved current fL. which is associated with 
the exact U(l) symmetry ( |3.10| ) we replace Xx exp{iWx)Xx and \ x — »■ x x exp(— iu x ) in the 



action ( |3.5|) and collect terms linear in tu x . But then the three link couplings in (|3.5| ) lead 
to an awkward looking divergence equation. The usual form of the divergence equation is 
obtained in the canonical model (|3.6|) . We find, 

d X = > 

jfMX = l^fJ-x^fix^iXxXx+fi + Xx+fiXx) ~ S U x£fj,uVtJ.X2(XxXx+fi ~ Xx+fiXx) ■ (4-9) 

The potential problem of the appearance of additional global symmetries on the lat- 



tice has been emphasized in particular in ref. [14|, in the context of fermion number non- 
conservation in the Standard Model. We shall assume here the following resolution, which is, 
as we believe, in accordance with current lore [15) . We can construct many currents, each of 
which reduces in the scaling region to a linear combination of the gauge invariant and the 
gauge variant e fJ-u A u . In particular, the exactly conserved but gauge non-invariant current 
j^ x will reduce to the divergence free combination jY — (i/ir)e^ u A u . Because this current is 
not gauge invariant, the corresponding conserved charge is unphysical. It may have a phys- 
ical (gauge invariant) component, but there is no reason why this should be conserved (for 
an exposition of the physics of the equivalent Schwinger model, see for example ref. ||16|| ). 

In this way the non-gauge invariance of our lattice model provides presumably a possible 
way out of the embarrassing exact global U(l) invariance. Another possibility to avoid 
difficulties with undesired global invariances is to construct the lattice models such that 



additional extra symmetries do not emerge j|, [TJ], |T7[. The potential problem of having a 
larger global symmetry group on the lattice than in the continuum target model, requires 
further detailed investigation. However, in this paper we shall restrict ourselves to a study 
of the currents J^ x and J A as given in the eqs. ( |4.6|) and ( |4.7| ). 

With the above definitions of the currents and using the staggered fermion formalism 
outlined in ref. [JlOj, we derived the lattice analogues of the Ward identities in sect. 2. 



Let us first concentrate on the vector current. The diagram in fig. |T]a gives a non-zero 
contribution, whereas the contribution from the typical lattice tadpole diagram in fig. |T]b 
happens to vanish. The amplitude then reads 

rVA/^ _ ,v*(p„-Pm)/2 f + tv«*s(q a )s(q v + p v ) + e ua s(q a + p a )s(q^ 

V [P) " h D(q)D(q + p) 

x c(^+^/2)c(g„ +p v /2)I[c(p j /2)c(q j +p j /2) , (4.10) 

j 

where s(qj) = sin(^), c(qj) = cos(^), D{q) = J2 a sin 2 q a + m 2 and J q = f^j^ d 2 q/n 2 . To 
calculate the continuum limit of (|4.10|) we let p and m approach zero, and separate the 
integration region into a small ball around the origin, |g| < 5, with radius 5 n/2 and the 
outer region, \q\ > 5 (see, for example, flOl ). We let <5 — > 0, m/8 — > and p/5 — > 0, with p/m 



fixed. In the inner region we can replace the integrand by its covariant form, while in the 
outer region the integrand is expanded in powers of m, p and only the non-vanishing terms 
are kept. For the integral (|4.10| ) the contribution from the outer region vanishes and only the 



contribution from the inner region remains, which is exactly a continuum loop integral with 
a spherical cut-off 6, and coincides with eq. ([2.5|), except for a factor of two corresponding to 
the two flavors. Consequently, after normalizing to one flavor we have the same contact term 
( pT8|) as in the cut-off regulated continuum theory. So we find the following vector current 
divergence relation on the lattice 

d;<JD = 2iF m /2K + 0(a), (4.11) 



where the contact term in the definition ( [4.6| ) is given in eq. fl2.8| ) and F x is a suitable 
form for the field strength, c.f. eq. ( |5.2|) below. Here d'^f^ x = Y^^(ffix~fx-p) is the divergence 
on the lattice and denotes the integration over the x fields and includes here and in the 
following also a normalization to one staggered flavor. The 0(a) indicates terms which arise 
due to the discretization and vanish when a —>■ 0. 

For the amplitude T^f we obtained the expression, 

Ttfip) = TX {a \v) + TX (b \p) + C£ , (4.12) 
J(p,~p^/2 [ m2( W + s(qjs(q u - p u ) ~ e^eypsjqcjsiqp + p p ) 

D(q)D(q + p) 

c(q»+pj2)c(q v +p v /2)\[c(q 3 +p 3 /2) 2 , (4.13) 



' J a 



(4 - 14) 

where T^f^ (p) and T^f^ (p) denote the contributions from the diagrams in fig. |]a and b. 
The amplitude of the tadpole diagram is independent of the external momentum p. After 
a normalization to one flavor, we find the lattice divergence relation for the axial-vector 
current, 

^(^) x = 2m(jf) x + 0(a), (4.15) 
where the contact term coefficient in (14.71) has to be chosen as 



1 



C = tS^ , r = - y— +h + I 2 j w 0.0625 , (4.16) 

j _ i r H^) 2 - s{ qi f)c{ qi Yc{ q2 f ^ Q2g3 _ i r S ( qi y _ 1 

1 2J q [ S ( gi )2 + s(g2 ) 2] 2 - ' 2 2 J q [ S ( qi y + s(q 2 ff 4" 

(4.17) 

The first two terms, 1/2-7T and I\, in (|4.16|) come from inner and outer region parts of the 
integral ( 4.13j ) and I 2 from the lattice integral (|4.14|) . 



An alternative form of the staggered fermion action is given by eq. ( |3.5| ), but with non- 
compact gauge fields, i. e. with the replacements c^ x 1 and s^ x A^ x . Of course, also 
this form reduces in classical continuum limit to the target model with two flavors. Using 



this modified form of the action we would have to drop the U fields from eq. ( |4.4| ). Then 
the tadpole diagram in fig. |l]b would not have given a contribution and we would instead 
have obtained C^f ~ —0.1875 5^ v which is larger than the result in ( |4.16|) . This shows that 



the form ( |3.5|) is more appropriate for restoring gauge invariance. However, unlike with the 
continuum point split current , a contact term is still needed because the staggered theory 
is not gauge invariant. 



5 Numerical results 



In this section we numerically compute the current divergence relations ([4.11 ) and ( 4.15| ) for 



smooth external gauge field configurations with variable amplitude. A similar test |18[ has 
been recently applied to a proposal using domain wall fermions [|19| . Since gauge invariance 
is violated in the staggered model it is interesting to investigate the effect of fluctuating 
gauge transformations on the divergence relations. As in an earlier work [^Oj we shall also 
investigate the effects of external gauge fields with Q ^ 1 . 

5.1 Smooth external gauge field configurations 



To test the relation ( [4.11| ) for the vector current, we use the same external fields as in ref. |JT8 



which are spatially constant, 

A lx = A sin(27rf/T) , A 2x = (5.1) 

with t = x 2 . We shall use here and in the following a lattice, with extents T and L in the 
time and spatial directions. The gauge fields are close to one everywhere provided the 
amplitude Aq is sufficiently small. The topological charge Q is equal to zero for this class of 
configurations. 

To reduce the discretization error, which arises naturally when transcribing a continuum 
gauge field configuration to the lattice, we use for F x in ( (4.11|) the average over the four 
lattice plaquettes adjoined to the point x: 

Fx = \Y, F i2*-b ( 5 - 2 ) 
4 b 

with F 12x = diA 2x - d 2 Ai x the plaquette field strength (d^f x = f x+fl - f x ). Fig. |a shows 
the divergence — id'^(J^ x ) ^ as a function of t for various values of A Q . The quantity F x /tt 
is represented in this plot by the full lines which were obtained by connecting the points at 
t — 1, . . . , T by a cubic spline. We use a lattice with T = L = 64 and periodic boundary 
conditions for the \ fields. To regulate the near zero mode we used a small non-zero mass 
m = 0.01. In fig. and in the following graphs we have multiplied the currents and the 
anomaly by the lattice volume L x T. For the smaller amplitudes the agreement between 
the anomaly and the current divergence is almost perfect, showing that the 0(a) effects in 
eq. ( |4.11| ) are very small. We find the maximal relative error to increase from 0.1% to 5% 
when Aq is raised from 0.1 to 0.25. 

To test the divergence relation ( |4.15| ) for the axial-vector current we use the field 



A lx = , A 2x = A sin(27rt/T) , (5.3) 

which in contrast to the previous case is longitudinal, i.e. e^ v d^A u = and d'^A^ ^ 0. In 
fig. ||b, ^{J^)^ is represented by the various symbols, which as in fig. ||a correspond to 

different values of the amplitude A . The quantity 2m(J x ) x is represented by the solid lines 
which were obtained by connecting the data points at t — 1, . . . , T by cubic splines. The plot 
shows that the relation fl4.15|) holds nicely within the given range of Aq values. The relative 



error is smaller than 6%. The complicated looking y4 -dependence of ^(J^) shows that 
for m/fl the effective action S e ff(A), which results after carrying out the x integration in 
the path integral, and thereto also (J^ x ) = —5S e ff(A)/SA lJiX , are non-linear functional of 
the vector potential A^ x . 



Figure 2: a) The divergence —id^(J^ x )^ as a function oft for several values of Aq and 

m = 0.01. The anomaly F x /n is given by the full lines, b) d'^(J^ x )^ as a function oft for 

several values of A and m = 0.01. The numerical results for 2m(J x ) x are represented by 
the full lines which were obtained by connecting the data points at various t by cubic splines. 



5.2 Fluctuating gauge degrees of freedom 



Since the model lacks gauge invariance it is instructive to investigate how relation Q4.ll ) 
changes if we perform gauge transformations 

-> n x u^n* x+ix , (5.4) 

on a smooth link configuration. To investigate the effect of Q field fluctuations, we can again 
compute the current divergence, but now averaged over an fl field ensemble, 

(W a = DSlfyjy)/™ , Z = J DQe^ , (5.5) 

with Boltzmann weight expS(fi). We will generate here the configurations of the gauge 
degrees of freedom in two different ways, depending on which formulation of the full model 
with dynamical gauge fields we are aiming for. 

As we mentioned in sect. 1, one possibility is to use non-perturbative gauge fixing and 
add counterterms to restore gauge invariance ||. With non-compact U(l) gauge fields (but 
still coupled compactly to the fermions), we may use a gauge fixing action — (/2 ^{^A^} 2 . 
The corresponding Fadeev- Popov determinant is independent of A^. Since A^ transforms 
into A^ x + d,j,uj x under a gauge transformation, this suggests to use the gauge fixing action 

S(Q) = {Ou x } 2 , Q x = expiu x , (5.6) 

X 

to generate the gauge degrees of freedom. Here □ denotes the lattice laplacian operator 
d'^d^ and ( is the gauge fixing parameter (£ = oo corresponds to the Landau gauge). The 
non-compact phases uj x are coupled through Q x = expiu x and the replacement (|5.4|) to 
the fermions. In fig. |3] the crosses represent the results for —i(d' fl J^ x ) n after averaging over 
1400 independent Q configurations at C — 10- The numerical result lies slightly below the 
solid line which represents here again the anomaly F x /ir for the given external gauge field 
configuration ( |5.1| ). This shows that the anomaly relation remains valid after multiplying the 
vector current by a factor which is slightly larger than one. Such a current renormalization 
is expected because there is no protection by symmetry. This result shows that a further 
numerical investigation of the gauge fixing approach may be technically feasible, at least for 
U(l) with non-compact gauge potentials where one has not to worry about the Fadeev- Popov 
factor. 

We find that at smaller values of £ the statistical fluctuations increase tremendously. 
Recall, however, that in the convention we are using the gauge coupling g is absorbed in 
A^ x , implying that ( oc 1/g 2 ; e.g. ( = 1/g 2 corresponds to the Feynman gauge. Since g 2 — ► 
in lattice units, large £'s are not unnatural. 

An alternative approach to regain a gauge invariant quantum model aims at a dynamical 
restoration of gauge invariance, by integration over all gauge transformations @, 0. In this 
approach the expected mass counterterm for the gauge field 2k Y^ X /i R- e U^ X) which transforms 
into 2nJ2x l j,R' e (^xU/xxil*^), suggests to use the action 

Xfl 

to generate the Q configurations in ( p.5| ). This action is identical to the action for the 
XY-model in two dimensions with hopping parameter k. The most ambitious scenario 0. || 
corresponds to choosing k in the vortex phase of the XY-model (in ref. denoted as 'scenario 



Figure 3: The divergence —i(d'^J^ x ) n as a function oft with the gauge degrees of freedom 
Q generated by the gauge fixing action fl5l% ) at ( = 10 (crosses) and by the XY-model action 
( |5. 7| j at k — 1.0 (circles). The anomaly F x /ir is represented by the full line. The dashed 
curve is obtained by a fit. 

C). In other scenarios the models would have to allow an interpretation of Q x as a (radially 
frozen) Higgs field, for which k would have to be close to the Kosterlitz-Thouless phase 
transition at k c « 0.5. 

For large values of k, deep in the spin wave phase, the fluctuations of Q are small and we 
expect that —i(d > IJt J^ x ) n = F x /tt within statistical errors. This is confirmed by a simulation 
at k — 2. For smaller values of n, when approaching the phase transition at k c « 0.5, 
the fluctuations increase dramatically. The result at k — 1.0 (indicated by circles in fig. |3D 
shows that the values for —i(d'^J^ x ) n now lie significantly below F x /tt (solid line). This 
result was obtained after averaging over 6400 independent Q configurations. The dashed 
line was obtained by fitting the numerical data to the ansatz — c(A 2ir/T) cos(2nt/T) with 
free parameter c. The good quality of the fit shows that the relation ( |4.11| ) remains valid after 
renormalizing the vector current by a factor 1/c ~ 1.20. In the vortex phase the fluctuations 
of —id' J^ x were so strong that even after an excessive increase of the statistics we were not 
able to get a meaningful estimate of the average current divergence. The danger is that 
the nice scaling behavior of the staggered fermions, which we could demonstrate for smooth 
external fields, is washed out when the gauge degrees of freedom fluctuate too strongly. 

An interesting difference between the actions (|5.6|) and ( |5.7|) is that in the latter case the 
small fluctuations in the spin wave phase are weighted by exp(KujOuj) , whereas in the former 
case a D 2 appears instead of the □, which leads to much smoother uj configurations. 



5.3 Configurations with Q ^ 



One would like to promote the divergence equation ( |2.4| ) to gauge fields with arbitrary 
topological charge. At first sight this leads to an apparent contradiction: On a torus with 
periodic boundary conditions for gauge invariant quantities, integration over the left-hand 
side of eq. ( |2.4|) gives zero, while the right-hand side gives 2iQ ^ 0. This paradox makes 
it interesting to see if the local divergence equation ( |2.4|) is valid also for topologically non- 
trivial gauge fields. 

In the massive Schwinger model the divergence equation for the axial- vector current reads, 
& JA, = 2mJ x + 2iq. After integration over the torus with periodic boundary conditions we 



fj, /IX 



find Q = mTr (75 [^)+ m) ) and the Atiyah-Singer index theorem is effectively valid p0 
A paradox is therefore avoided by the fermion mass term in the massive Schwinger model. 
The paradox in the axial- vector model would not emerge if we would replace the usual mass 
term in eq. ( |2.1| ) by a Majorana mass term which makes the model equivalent to the massive 
Schwinger model. 

However, it is not as easy to deal with Q 7^ in the axial-vector model as previously 
in the Schwinger model (see ref. [Q). The reason is that we recover gauge invariance only 
if the gauge potentials are smooth (and, of course, m = 0). For example, the perturbative 
derivation of the divergence relations is valid only when the momenta of the external gauge 
field are negligible compared to the cut-off. On a periodic lattice a typical lattice gauge 
field with non-zero topological charge is not smooth. As an example we investigate the effect 



of configurations with constant F\2 X used in the vector case in ref. [20] 



We consider a configuration with a constant field strength F = Fi2 X = 2ttQ/TL, 

U lx = exp(iFt) , t = 1, . . . ,T , 

U 2x = l, t = l,...,T-l, (5.8) 

U 2x = exp(iFTxi) , t = T . 

For small Q the link field is close to one and smooth everywhere except for t = T where 
U2 contains a transition function [pCfl . When shifting t from T — 1 to T and then to T + 1 = 1 
(mod T) U2 makes a jump as a function of t (e.g. for Q = 1 and x\ — L/2 U 2 jumps from 
+ 1 to —1 and then to +1 again). Also U\ jumps when shifting t from Ttot = T + l = l 
(mod T). In gauge invariant models these transition functions are invisible. 

We expect the divergence equation to hold with small 0(a) corrections in the region of 
the lattice where the gauge fields are sufficiently smooth, but expect deviations in the region 
near the transition function. As in fig. |2] we have plotted in fig. [| the divergence —id'^(J^ x ) 
(squares) for the space slice x\ = L/2 as a function of t/T with L = T = 64 (diamonds) 
and L = T = 32 (crosses), but now for the gauge field configuration (|5.8|). We used here 
antiperiodic boundary conditions for the x fields and m = 0. For other space slices we got 
similar plots. The solid line represents again F x /tc, which now is independent of x. Far away 
from the time slice which carries the transition function the agreement of — id , fl (J'^.) with 
F x /tt is satisfactory, however in the vicinity of this time slice the deviations become huge 
(the values at t = 31, 32 (T = 32) and t = 62, 63, 64 (T = 64) were dropped from the graph 
since they are much larger than 20). The figure indicates also that the region of disturbance 
shrinks for increasing lattice size. The strong deviations induced by the transition function 
is analogous to the enormous fluctuations we observed in —i&^J^) when lowering k in 
the action ( |5.7| ). The increase of the vorticity renders the effective gauge field configuration 
Q x U flx Ql + f i very rough, similar to the above gauge field at t — T. 

The above results for the Q = 1 configuration are not entirely satisfactory. To avoid the 
large errors with topological non-trivial gauge fields we will have to follow the mathematicians 



Figure 4: The divergence — id'^(J^ x )^ as a function oft/T for a configuration with Q = 1 
on two different lattices (T = L = 32, 64). We have used m = 0, and antiperiodic boundary 
conditions for the x field. The solid line represents the anomaly F x /n. 

and introduce charts in which the gauge potentials are smooth, such that the Dirac operator 
has negligible discretization errors in any open region (see e. g. ref. [pl|). A proper treatment 
of this will have to be done in the future. 

6 Discussion 

We have showed that the staggered fermion model ( |3.5| ) can reproduce continuum Ward 
identities after incorporating the appropriate counterterms. The coefficients of these coun- 
terterms have been computed in this paper within lattice perturbation theory. Using smooth 
external gauge fields with zero topological charge we have numerically verified the validity 
of the divergence relations for the vector and axial-vector currents and found good accuracy, 
for gauge field configurations with jL 7 ^ — 1|^0.2. 

We tested the sensitivity of the divergence relation for the vector current to fluctuating 
gauge degrees of freedom which were generated either with the gauge fixing action (|5.6|) or 
the scalar field action (|5.7|) . After renormalizing the currents with a finite factor > f, the 
divergence relation remains valid in both cases, provided that the fluctuations of the gauge 
modes Q x are not too strong. This is an encouraging result for the gauge fixing approach to 
chiral theories. When the gauge modes becomes less constrained, the induced fluctuations 
are very severe and might wash out the anomaly signal completely. 

The disconcerting effect of rough gauge transformations is also seen when the current 
divergence is measured for a non-smooth gauge field with topological charge one, used in 
earlier tests in QED2. We found the divergence relation for the vector current to be strongly 
violated near the region of the lattice where the gauge potentials lack smoothness. This 



means that in a description which violates gauge invariance at the cut-off level, topologically 
non-trivial gauge fields have presumably to be dealt with the full apparatus of charts and 
transition functions that are smooth in space and time. 

The results of this paper show that our staggered fermion approach passes a first test 
in reproducing the current divergence relations. It is not clear yet at this stage whether it 
is possible to obtain a valid quantum model after the integration over all gauge field con- 
figurations has been carried out in the path integral. This question shall be addressed in a 
separate publication JjJ. Further clarification is also needed of the question how (if indeed) 
the additional global U(l) invariance on the lattice does not give rise to a local conservation 
law. 
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